Search for: All records

Volcano monitoring and eruption forecasting require accurate characterization of transcrustal magmatic structures to place volcanic unrest in context within the system where it occurs. Structural imaging using local seismicity is limited to seismogenic depths. Here, we exploit arrivals in teleseismic receiver functions that change polarity with backazimuth to image two surfaces beneath Akutan volcano in the Aleutian arc. The two surfaces delineate an upper to midcrustal inverted conical volume that deepens and thickens away from the volcanic center, with thicknesses of 3–13 km. The top of the volume is at depths of 2–3 km below sea level at distances of ∼5–15 km from the caldera center. The bottom is at depths of 7–15 km at the same distances, and the cone’s thickness increases outward from ∼5 to ∼10 km. The signal is best fit by a volume with anisotropy with fast symmetry planes that dip outward from the center and downward increases in shear velocity at both interfaces. The upper boundary coincides with the top of Akutan’s volcanotectonic (VT) seismogenic zone, with the VT seismicity exhibiting outward dipping planar features that match the anisotropic fast plane orientation within the volume. The bottom of the anisotropic volume is below the termination depth of the majority of the VT seismicity and is therefore likely associated with the brittle–ductile transition. Long-period (LP) events associated previously with magma movement are concentrated below the anisotropic VT volume. Because of the strong spatial association with VT seismicity, we interpret the volume as consisting of concentric outward dipping faults and dikes that align the seismogenic response to stress changes from magmatic processes. Our observations map this volume independent of the present-day seismicity distribution and thus provide a spatially more complete image of the magmatic system. 

Abstract Megathrusts at convergent plate boundaries generate the largest and some of the most hazardous earthquakes on Earth. However, their physical properties, including those influencing fault slip accumulation and release and earthquake‐related surface displacements, are still poorly constrained at critical depths. Here, we combine seismic imaging and geodetic modeling to investigate the structure and mechanical behavior of the Main Himalayan Thrust fault (MHT) in the center of the 2015 Mw 7.8 Gorkha rupture in Nepal. Our results from two independent observations consistently suggest the presence of a channel associated with the MHT with high compliance (shear modulus as low as ∼4 GPa) and strain anisotropy (stiffer in the vertical orientation than in the horizontal), likely arising from a weak subducting layer with north‐dipping foliation. Such mechanical heterogeneity significantly influences the quantification of short‐term fault kinematics and associated earthquake potential, with implications on across‐scale dynamics of plate boundaries in Himalaya and elsewhere. 

Abstract A long‐standing question is how felsic continental crust is differentiated from its mafic parent mantle magmas. One currently proposed fundamental mechanism is lithospheric foundering and loss of dense material into the mantle. A type locality is the young extinct arc forming the Sierra Nevada, California. Here, we image a distinct anisotropic shear layer below the crust‐mantle boundary using receiver functions. The sense of shear is consistent with west‐ to southwestward removal of lithosphere. The shear signal is strongest in the southern Sierra, where lithospheric foundering was proposed to have concluded several million years ago, and is deeper and less pronounced in the central Sierra, where ongoing lithospheric foundering is corroborated by a band of unusually deep (40+ km) seismicity along the western foothills. Our observations provide progressive snapshots of a lithospheric foundering process spanning several million years and hundreds of kilometers, illuminating a fundamental differentiation process by which continents are built. 

The composition of the crust is one of the most uncertain and controversial components of continental estimates due to (1) limited direct access and (2) inconsistent indirect assessments. Here we show that by combining high-resolution shear velocity (Vs) models with newly measured with newly measured ratio of compressional wave velocity (Vp) and Vs, or Vp/Vs ratio, for the crystalline crust, a 3-D composition (SiO2 wt%) model of the continental crust can be derived with uncertainty estimates. Comparing the model with local xenolith data shows consistency at mid and lower crustal depths. The spatial patterns in bulk SiO2 content correlate with major geological provinces, including the footprints of Cenozoic and Mesozoic mafic volcanism in the western U.S., and offer new insight into the composition and evolution of the continental U.S. 

Abstract As North America collided with Africa to form Pangea during the Alleghanian orogeny, crystalline and sedimentary rocks in the southeastern United States were thrust forelandward along the Appalachian décollement. We examined Ps receiver functions to better constrain the kinematics of this prominent subsurface structure. From Southeastern Suture of the Appalachian Margin Experiment (SESAME) and other EarthScope stations on the Blue Ridge–Piedmont crystalline megathrust, we find large arrivals from a 5–10-km-deep converter. We argue that a strong contrast in dipping anisotropic foliation occurs at the subhorizontal Appalachian décollement, and propose that such a geometry may be typical for décollement structures. Conversion polarity flips can be explained by an east-dipping foliation, but this orientation is at odds with the overlying northeast-trending surface tectonic grain. We suggest that prior to late Alleghanian northwest-directed head-on collision, the Appalachian décollement accommodated early Alleghanian west-vergence, independent of the overlying Blue Ridge–Piedmont structural inheritance. The geophysical expression of dipping anisotropic foliation provides a powerful tool for investigating subsurface kinematics, especially where they are obscured by overlying fabric, to disentangle the tectonic complexities that embody oblique collisional orogens. 

Abstract Seismic anisotropy constitutes a useful tool for imaging the structure along the plate interface in subduction zones, but the seismic properties of mafic blueschists, a common rock type in subduction zones, remain poorly constrained. We applied the technique of electron backscatter diffraction (EBSD) based petrofabric analysis to calculate the seismic anisotropies of 14 naturally deformed mafic blueschists at dry, ambient conditions. The ductilely deformed blueschists were collected from terranes with inferred peak P‐T conditions applicable to subducting slabs at or near the plate interface in active subduction zones. Epidote blueschists display the greatestPwave anisotropy range (AVp ∼7%–20%), while lawsonite blueschist AVp ranges from ∼2% to 10%.Swave anisotropies generate shear wave splitting delay times up to ∼0.1 s over a thickness of 5 km. AVp magnitude increases with glaucophane abundance (from areal EBSD measurements), decreases with increasing epidote or lawsonite abundance, and is enhanced by glaucophane crystallographic preferred orientation (CPO) strength. Two‐phase rock recipe models provide further evidence of the primary role of glaucophane, epidote, and lawsonite in generating blueschist seismic anisotropy. The symmetry ofPwave velocity patterns reflects the deformation‐induced CPO type in glaucophane—an effect previously observed for hornblende on amphibolitePwave anisotropy. The distinctive seismic properties that distinguish blueschist from other subduction zone rock types and the strong correlation between anisotropy magnitude/symmetry and glaucophane CPO suggest that seismic anisotropy may be a useful tool in mapping the extent and deformation of blueschists along the interface, and the blueschist‐eclogite transition in active subduction zones. 

Abstract Azimuthal variations in receiver function conversions can image lithospheric structural contrasts and anisotropic fabrics that together compose tectonic grain. We apply this method to data from EarthScope Transportable Array in Alaska and additional stations across the northern Cordillera. The best-resolved quantities are the strike and depth of dipping fabric contrasts or interfaces. We find a strong geographic gradient in such anomalies, with large amplitudes extending inboard from the present-day subduction margin, the Aleutian arc, and an influence of flat-slab subduction of the Yakutat microplate north of the Denali fault. An east–west band across interior Alaska shows low-amplitude crustal anomalies. Anomaly amplitudes correlate with structural intensity (density of aligned geological elements), but are the highest in areas of strong Cenozoic deformation, raising the question of an influence of current stress state. Imaged subsurface strikes show alignment with surface structures. We see concentric strikes around arc volcanoes implying dipping magmatic structures and fabric into the middle crust. Regions with present-day weaker deformation show lower anomaly amplitudes but structurally aligned strikes, suggesting pre-Cenozoic fabrics may have been overprinted or otherwise modified. We observe general coherence of the signal across the brittle-plastic transition. Imaged crustal fabrics are aligned with major faults and shear zones, whereas intrafault blocks show imaged strikes both parallel to and at high angles to major block-bounding faults. High-angle strikes are subparallel to neotectonic deformation, seismicity, fault lineaments, and prominent metallogenic belts, possibly due to overprinting and/or co-evolution with fault-parallel fabrics. We suggest that the underlying tectonic grain in the northern Cordillera is broadly distributed rather than strongly localized. Receiver functions thus reveal key information about the nature and continuity of tectonic fabrics at depth and can provide unique insights into the deformation history and distribution of regional strain in complex orogenic belts. 

Abstract Seismic anisotropy is controlled by aligned rock‐forming minerals, which most studies attribute to solid‐state shear with less consideration for magmatic fabric in plutonic rocks (rigid‐body rotation of crystals in the presence of melt). Our study counters this traditional solid‐state bias by evaluating contributions from fossil magmatic fabric. We collected samples from various tectonic settings, identified mineral orientations with electron backscatter diffraction and neutron diffraction, and calculated their bulk rock elastic properties. Results indicate that magmatic fabric may lead to moderate to strong anisotropy (3%–9%), comparable to solid‐state deformation. Also, magmatically aligned feldspar may cause foliation‐perpendicular fast velocity, a unique orientation that contrasts with a fast foliation typical of solid‐state deformation. Therefore, magmatic fabric may be more relevant to seismic anisotropy than previously recognized. Accordingly, increased considerations of magmatic fabric in arcs, batholiths, and other tectonic settings can change and potentially improve the prediction, observation, and interpretation of crustal seismic anisotropy. 

Abstract Earthquakes are known to occur beneath southern Tibet at depths up to ∼95 km. Whether these earthquakes occur within the lower crust thickened in the Himalayan collision or in the mantle is a matter of current debate. Here we compare vertical travel paths expressed as delay times between S and P arrivals for local events to delay times of P-to-S conversions from the Moho in receiver functions. The method removes most of the uncertainty introduced in standard analysis from using velocity models for depth location and migration. We show that deep seismicity in southern Tibet is unequivocally located beneath the Moho in the mantle. Deep seismicity in continental lithosphere occurs under normally ductile conditions and has therefore garnered interest in whether its occurrence is due to particularly cold temperatures or whether other factors are causing embrittlement of ductile material. Eclogitization in the subducting Indian crust has been proposed as a cause for the deep seismicity in this area. Our observation of seismicity in the mantle, falling below rather than within the crustal layer with proposed eclogitization, requires revisiting this concept and favors other embrittlement mechanisms that operate within mantle material.